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ABSTRACT 

Single-armed, stationary density waves can exist even in disks with only weak self- 
gravity, provided that the rotation curve is dominated by a central mass. Such waves 
could play a significant role in the transport of angular momentum. By variational 
methods, we derive nonlinear versions of the dispersion relation, angular momentum 
flux, and propagation velocity in the tight-winding limit. The pitch angle increases 
with amplitude until the tight-winding approximation breaks down. By other methods, 
we find a series of nonlinear logarithmic spirals which is exact in the limit of small 
disk mass and which extends to large pitch angle. These waves may be supported by 
low-mass protoplanetary disks, and perhaps by compact molecular disks in galactic 
nuclei. 

Key words: accretion, accretion discs - galaxies: nuclei - hydrodynamics - Solar 
system: formation - waves 



1 INTRODUCTION 

The number of spiral arms in selfgravitating disks is influenced by the rotation curve. As observed by Lindblad, near-harmonic 
potentials and linearly rising rotation curves favor two-armed spirals: free-particle orbits are slowly precessing ellipses centered 
on the minimum of the potential; spiral wav es can be built up from such orbits with relatively weak interactions needed to 
force them to precess at a common rate (cf. Binney fc Tremaine (1987) ). Since galactic potentials are often approximately 
harmonic at small radii, this construction may in part explain the prevalence of two-armed spirals and bars in disk galaxies. 

The rotation curves of astrophysical disks smaller than galaxies are usually dominated by a central point-like mass: 
examples include accretion disks, protoplanetary disks, planetary rings, and masing molecular disks recently discovered in 
galactic nuclei (Section 4 and references therein). Free-particle orbits in such disks are slowly precessing keplerian ellipses. In 
this case, Lindblad's construction favors single-armed spirals, again with small angular pattern speed. Given the astrophysical 
importance of keplerian disks, it is unfortunate that single-armed spirals have received much less theoretical attention than 
double-armed ones. The visual fascination of beautiful spirals in galactic disks, and the absence of comparably well resolved 
images of keplerian disks, may partly explain the neglect. 

The present paper is therefore dedicated to the propagation of single-armed spirals in nearly keplerian gaseous disks. We 
concentrate on the limit where the self-interaction of the disk is weak compared to its interaction with the central mass. This 
implies small or vanishing angular pattern speeds. The nonlinear regime is emphasized because it is gratifyingly tractable, 
and more importantly because one would like to know the maximum angular momentum flux that can be transmitted by 

Shu 1989| ; fchu et al. 199C| ; 



density waves. We do not study the linear instabilities of single-armed waves (cf. 



Heemskcrk, Papaloizou, & Savonije 1992) or their excitation by companion masses (Yuan 



Adams, Ruden, 

Cassen 198^ )7 



In Section 2, we adopt the familiar tight-winding approximation and extend it to nonlinear spirals. This has already 
been done using other methods, notably by Shu, Yuan, fc Lissauer (1985)| ; 5hu, et al. (1985)| ; and Borderies, Goldreich, & 
Tremaine (1985, 1986). Our mathematical approach, however, is adapted from the variational methods of Whitham (1974), 



which are particularly efficient for describing the propagation of a nonlinear wave train outside the region where it is excited, 
especially if dissipation is weak. In addition to the nonlinear dispersion relation and angular momentum flux, the variational 
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approach provides the nonlinear analog of the group velocity, which as far as we know have not otherwise been derived for 
density waves (Appendix B). Also, previous work has mainly been concerned with nonlinear extensions of the long branch of 
the WKB dispersion relation, whereas we are primarily concerned with the short branch. 

We show in Section 2 that for single-armed waves in keplerian potentials, self-gravity is best measured not by the familiar 
Q parameter (Toomre 1964), but by the ratio (^) of "Jeans length" (Aj) to disk radius. In fact, almost all of our analyses 
pertain to the limit as Q — > oo at finite Xj/r. 

The results of Section 2 show that single-armed waves unwind with increasing nonlinearity. Therefore Section 3 presents 
nonlinear logarithmic spirals that do not depend upon the tight-winding approximation. The method of their construction 
assumes radial self-similarity, which restricts us to disks whose azimuthally averaged surface density scales with radius as 
j- -3 / 2 . "We believe, however, that many of the results are more generally valid. Solutions do not seem to exist unless Xj/r <> 4. 

We summarize our results in Section 4 and discuss applications to astrophysical disks. In particular, we show that the 
minimum-mass solar nebula was probably in the regime considered here, namely large Q but moderate Xj/r, so that stationary 
single-armed waves may have been important. 



2 TIGHTLY- WRAPPED WAVES 

In this section, we consider single-armed waves of small pitch angle. In Section 2.1, we summarize the linear theory, emphasizing 
those aspects that are peculiar to m = 1 waves in nearly keplerian potentials. These results provide background and guidance 
for the development of the nonlinear theory in Section 2.2, which itself has two parts: first we sketch a variational formalism 
for nonlinear wave trains in general, and then we apply this formalism to tightly-wrapped spiral density waves. 



2.1 Linear Theory 

The dispersion relation for tightly wrapped m = 1 spiral waves is 

(cj-n j 2 = k 2 -2ivGEk + c 2 k 2 . (1) 

However it is excited, as the wave propagates inward, the pattern speed (ui) becomes negligible when compared to the local 
angular velocity of the disk (f2). We study the stationary limit, uo <C Q, which is appropriate to waves far inside their corotation 
radius close to a central point mass. In this region, both a 2 and Q, 2 are dominated by contributions from the central mass, 
but these large terms cancel one another to leading order in equation (|l|), leaving a residual due only to the disk. 

To simplify the analysis we consider a razor thin disk and take a polytropic equation of state P — KTP . For a logarithmic 
spiral pattern, the density and temperature profiles should be power laws, £ oc r - ' 3-1 and c 2 oc r _/3 where < f3 < 1. This 
requires 7 = (1 + 2/3)/(l + /3), so 7 can range from 1 to 3/2. Now k 2 — Q. 2 , the self-gravity term, and the pressure term 
in the dispersion relation all have the same scaling with r if the pitch angle cot _1 (fcr) ~ 1/kr is constant. Introducing the 
dimensionless parameters 



(2) 



2_ C 2 _X^ 

a 2ivGZr 2nr ' 
and 

g(p) = ( * ~ " 2)r ^ 2 +0.7536/3(1 -/3)- 2a 2 (l-/3 2 ) if < < 1. (3) 

7TG-Zj 

we write the dispersion relation in dimensionless form 

g - 2\k\r + 2a 2 (kr) 2 = 0. (4) 

Notice that k 2 — Q 2 is independent of the central point mass. The term involving a 2 in eq. (^) expresses the influence of pressure 
on the disk rotation curve; the rest is due to the disk's self-gravity. For small values of a 2 , the pitch angle cot _1 (fcr) « a 2 for 
the physically relevant root of the dispersion relation (see Section 3.1), so the re 2 — fl 2 term is negligible. 

To see how the amplitude of the linear wave varies with radius we use the fact that wave angular momentum is conserved. 
The angular momentum flux due to spiral waves arises from g ravitational and pressure tor ques, and from advective transport. 
The advective transport is proportional to (1 — c 2 |fe|/27rGE) ( |Goldreich fc Tremaine 1979 ) which vanishes to leading order in 



1/kr for these waves. The flux due to gravitational and pressure torques can be calculated for tightly- wound patterns in the 
WKBJ approximation: 

r<5<& 2 7r 2 rG5E 2 

Ywkb = sign(fc)^- = sign(fc) — (5) 
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Hence dT WKB /dr = implies <5E/E oc r' 3- ' 1 '' 2 ', so that the fractional amplitude of the wave increases inward if j3 < 1/2. 

The m — 1 waves have another special feature worth remarking upon here: in principle, they may exert a net force on 
the central mass. The component of the force along direction 8 in polar coordinates is 

27T OO 

G<5£(r, 0') 



,■2 



■cos(0 -6')rdrd6. 



o o 



If <5E(r, 9') oc exp(im9'), then the integral above can be nonzero for m = 1. In all cases of interest to us, however, 5E will be 
an oscillatory (wavelike) function of radius, so that the force is dominated by the endpoints of the radial integration. Thus, 
the displacement of the central mass cannot be included in a local analysis such as ours, whether linear or nonlinear, and we 



shall neglect it. Fo r a careful global treatment of this issue in nonlinear as well as linear regimes, see Heemskerk, Papaloizou 



Sa^ onije (1992) . These authors find that modes capable of displacing the central mass tend to be evanescent in radius 



rather than wavelike. 



2.2 Nonlinear Theory 

Approximate solutions for a wide variety of nonlinear wave trains arising from mathematical physics can be found with 
remarkable efficiency by varational methods. Since these methods are not widely used in astrophysics, we digress briefly to 



summarize the main ideas. For a more complete exposition, see Whitham (1974) . 
The variational approach rests on the following assumptions: 

(A.l) The fundamental equations of motion are derivable from an action principle, hence generally nondissipative. 
(A. 2) The wavelength and wave period of interest are small compared to the length scale and time scale of variations in the 
background. 

(A. 3) The waves are locally periodic in time and space, at least approximately. 

Assumption (A.l) can be relaxed by tacking dissipative terms onto the equations of motion at a late stage in the 
procedure, but we will not do so here. Assumptions (A. 2) and (A. 3) are shared with linear WKBJ theory, but since the 
principle of superposition does not hold, one cannot synthesize a solitary wave — a highly localized wave packet — by linearly 
combining periodic wave trains. On the other hand, the third assumption permits gradual changes in the wavelength, wave 
period, and wave amplitude, provided that they occur on scales long compared to the wavelength and period. 

With these assumptions, the following procedure accurately approximates nonlinear wave trains. The dynamical variable 
(it) depends in the first instance on one or more independent spatial variables (xi, . . . ,x n ) and on time (t); u might in fact 
represent multiple dependent variables. By assumption (A.l), u obeys an action principle: 

8 I L{u,dtu,d\u, . . . ,d n u;xi, . . . ,x n ,t)J(x)d n xdt — 0. (6) 



Here J(x) is just the Jacobian of the transformation from (xi, . . . ,x n ) to cartesian coordinates. We have indicated that L 
depends on first spatial derivatives of u, but in fact it may depend on spatial derivatives of any order, as is the case for 
self-gravitating density waves. For solutions of the action principle obeying assumption (A. 3), u can be written as a periodic 
function of a single phase (*!'): 

u = £/(*) = [/(* + 2tt), 9 = 9(xi,...,x n ,t). (7) 
The frequency and wavenumber are defined as derivatives of 9: 

cu = -dt^, fesWs— . (8) 

Hence d t u — > —uoU' and diU — > kiU 1 . The critical approximation is now to average over the phase 9 while pretending that 
gj, ki, , k„ and t, X\ , . . . x„ are constant: 

L(U,U',u,ki,...;xi,...,t) = J L{U,-ojU',kxU',...;xx,...,t)J{x)^. (9) 

o 

The action principle (^) reduces to 

8 / L{U,U';iu,k 1 ,...-x 1 ,...,t)J(x)d n xdt = 0, (10) 
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where S now indicates variation of U(ty) and ^(xi, . . . ,x n ,t). 

Stationarity of the action with respect to ty, which enters only through its derivatives 



d_ dL 
dt \ duo 



5^ dxt y dh 

8=1 X 



0. 



implies the conservation law 



(11) 



The quantity OL/dco is the density of wave action, and ~dL/dki is the i component of the corresponding flux. 

Stationarity of the action with respect to U yields both the nonlinear dispersion relation and also equations for the 
functional form of [/(*) [cf. [Whitham 1974 . 

To apply this formalism to spiral density waves, the radial wavelength must be small compared to the radius in the disk 
(r), since the surface density and other background properties vary significantly over radial distances comparable to r. For 
single-armed spirals, a short radial wavelength implies that the wave is tightly wrapped. So we may as well take advantage of 
tight winding to simplify our formulae. 

For dynamical variables, we take the radial and angular displacements of a fluid particle from a circular orbit: 



5r = r p (t), 60 = p {t)-n(a)t-0 p (0), 



(12) 



where [r p (t) , 8 p (t)] are the polar coordinates of the particle at time t. Here a and 8 P (0) are labels that follow the particle, 
but we will soon interpret a as the semimajor axis of a keplerian orbit that closely approximates the particle's trajectory. 
Q(a) is the mean angular velocity of that trajectory, which is not exactly the same as the mean motion of a keplerian orbit 
of semimajor axis a, because the pressure and self-gravity of the disk modify its rotation curve. For independent variables we 
take (a, 0, i). In tightly- wound waves where the streamlines do not actually cross, 5r/a 1 but dSr/da may approach unity. 
Therefore in constructing the lagrangian, we keep only the lowest important order in Srja but all orders in dSr/da. 

In the absence of any collective effects, the lagrangian and the action would be those of a collection of free particles: 



I/free (Sr, Sr , a) adadOdt, 



fdSry 

\~dTJ 



k (Sr) 



So (a), 



(13) 



where 



2a 3 da ^ > ^ 
is the epicyclic frequency, and So (a) is the surface density that fluid elements would have in the unperturbed axisymmetric 
disk. The corresponding Euler-Lagrange equation for each fluid element is Sr = — n 2 r. 

To represent the self-interaction of the disk, we subtract from 7f ree the time integral of two potential-energy terms. The 
first of these is the thermodynamic internal energy, VKmt- Absent dissipation, the two-dimensional pressure and density of a 
given fluid element scale as P oc E 7 , and the internal energy per unit mass is jP/(j — 1)E. Let Po(a) be the 2D pressure that 
a fluid element would have in a circular orbit. Its actual surface density is 



1 + 



dSry 
~da~) 



■£(a,8,t) = Eo(a) (l 

in the tight-winding approximation. It follows that the total internal energy of the gas is 
eg (a) 



Hi 



7(7 - 1) 



/ dSr\ 
+ 



"(7-1) 



Ho(a)adad9, 



(15) 



(16) 



where Cq(o) = 7P0/E0 is the square of the unperturbed sound speed. 

The gravitational self-energy of a tightly-wound wave is accurately expressed by the "thin-wire" approximation. On the 
scale of the radial wavelength, curves of constant surface density are approximately straight lines. The gravitational interaction 
energy of two lines of mass per unit length (/Ui,/i2) and separation S12 is Gfii/J.2 In | S12 1 per unit length. We put hi — Eo(a)dai, 
/12 = Eo(a)da2, where a = (a± + «2)/2, since Eo(a) varies only on the scale of a. And since the pitch angle of the streamlines 
is small, S12 = \a\ + Sri — (a.2 + <5?*2)|, where Sri = r(cn, 8, i) and 6 = (9i + #2)/2. Therefore, 



G d8 aT, (a)da / d(a\ — 02) In 



Sri — Sr-2 



ai — ct2 



(17) 



We have replaced S12 by Si2/(ai — 02) inside the logarithm; this changes W gra , v by a term independent of Sr, which has no 
effect on the equations of motion. The integration over (ai — 02) is also insensitive to its limits, since the dominant contribution 
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comes from |oi — ai\ < a in a tightly-wound wave. The above expression for Wgrav describes local self-gravity but not the 
long-range contribution of the unperturbed axisymmetric disk surface density to the rotation curve. The latter is included in 
the mean angular velocity Q(a) and in the epicyclic frequency (jl4j) derived from it. 

For single-armed wavelike solutions, Sr and 59 depend on (a, 9, t) via the phase variable 



* = 6- 



- ojt. 



In principle the wavenumber has two components, but since the angular component kg 

t _ a* , 

k = — = -<j> (a) 
Oa 



(18) 

8^/89 is always unity, we write 

(19) 



for the radial component. 

Our next step is to average the full lagrangian over one period in <]/. This task is simplified by the fact that we are 
interested in situations where the interior mass of the disk is small compared to that of the central object (M) . We introduce 
a small parameter 



2na So (a) 



(20) 



where a is typical of the radii of interest. The self-interaction terms Wi n t and W gla , v are formally of order r\ compared to the 
free action (|l3|). To leading order in 77, the trajectory of a fluid element must therefore be a solution of the Euler-Lagrange 
equation derived from L{ lac alone. Such a solutions is sinusoidal in t, and since it must depend on t through ty, it can be 
written in the form 



Sr(a, 9) = ae(a) cos ty, 



(21) 



where e is a dimensionless measure of the amplitude. For r\ = 0, the only forces would be those of an exactly keplerian potential; 
in that case, ui = 0, and e and <j> could be chosen as arbitrary constants independently for each fluid element. For r\ 7^ but 
small, e and <f> must be smooth functions of a because the self-interaction terms would diverge if neighboring streamlines were 
to intersect. By substituting the form ( |2l| ) into the lagrangian and applying the phase-averaged action principle, one finds 
relations between e, u, and k to leading order in r\. At higher orders in 77, there are corrections to these relations and to the 
functional form (^l[). 

Note that the small parameter r\ is independent of the degree of nonlinearity of the wave. The latter is measured in the 
tight-winding approximation by the streamline-crossing parameter 



qo ~ \k\ae 



(22) 



Streamlines intersect if and only if go > 1, as can be seen by differentiating eq. (^]) with respect to a assuming dine /da — 
0(1/ a) <Si k, and then referring to eq. (|l5|). We will assume that qo < 1, but not necessarily go -C 1, since otherwise dissipation 
is inevitable. Hence since ka 3> 1 is necessary for the tightwinding approximation, we must have e< 1. 
The time derivative in eq. ( ]l3| ) is the total time derivative following a fluid element, 

dSr dSr dSr 

— ^ — + n {a )—. 

The angular velocity 9 = Q + 0(e), but we ignore the O(e) correction when it multiplies angular derivatives, since these are 
small compared to radial derivatives. Hence the phase average of Lf ree is 



Z7T 

f d* E (a) r 2 2] 2/ n 

/ ° 27r = — 4 — ^ ~ ~ K -I 6 



(23) 



The phase averages of the internal and gravitational energies are derived in Appendix A, and are expressed in terms of 
the following integral functions: 



U{qo) 



and 



W(qo) 



7(7- 



7r|go| 



dy In 




1 - 



sin(goy) 



(24) 



(25) 
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Figure 1. Nonlinear dispersion relation for stationary single-armed waves at a 2 = 0.2, /3 = 1/2, and 7 = 4/3. Abscissa: streamline- 
crossing parameter qg (eq. [p2|) or q (eq. pc[]). Ordinate: dimensionless radial wavenumber ka or a, equal to cotangent of pitch angle. 
Solid curve: Short branch in tight- winding theory (ka vs. qo). Dashed curve: Unphysical long branch. Filled Points: Logarithmic spirals 
(a vs. q). Hollow Point: Linear theory ( Section 3.1) 



These functions are normalized so U — 1 + 0(g 2 ) and W = 1 + 0(q^). Both increase monotonically with go and both are 
finite but nondifferentiable at go = 1: W '(go) oc (1 — go) 1 ^ 2-7 as go —* 1, while W'(go) oc ln(l — go). 
The phase-averaged action is therefore 



L(e,w,k,a)2Tradadt. (26) 
It is convenient to divide L into a part Lo that is independent of u and a residual: 

L(e,u,k,a) — Lo(e,k,a) + — j— ^ w[w — 2S7(a)]a 2 e 2 , 
7rGEoa 



Lo(e, fc, a) 



ige 2 - |go|W(g )e + o- 2 g, 2 W(g ) 



2 

The averaged action (fed) must be stationary with respect to e; imposing dL/de — yields 



(27) 



5 e 2 - C sign(fc)-^-[g 2 W(g )] + a 2 g ^[g 2 W(go)] - w(w - 20)-^=- = 0. (28) 
ago ago 7rG2jo 

In the limit e 1 at fixed ka, every term in this expression is proportional to e 2 , and we recover the linear dispersion relation 
(@) or (@)j since W(0) = W(0) = 1 and W(0) = W'(0) = 1. In general, however, the dispersion relation relates the amplitude 
of the wave to its wavelength and frequency. If we fix one of these three quantities, then the other two are constrained to a 
curve. In this paper, we are most interested in u> = 0. Figure [j] displays some representative dispersion curves for stationary 
waves. We have used go rather than e as a measure of wave amplitude because it more directly controls the importance of 
pressure and gravitational forces. 

The Euler-Lagrange equation for $ is (|j"l|), with the jacobian J — 2ira in this case. Hence the wave action density (action 
per unit physical area) is 

_ 81 SpaV 

P =d^ = {LJ - Q) ^2— (29) 
This can also be interpreted as the angular momentum density of the wave. The density is negative for stationary waves 
(u> — > 0). The corresponding radial flux is 
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Figure 2. Wave torque (azimuthally integrated radial angular momentum flux) for cr 2 = 0.2, fi = 1/2, 7 = 4/3. Solid curve: Tight- 
winding theory; V = 2-naf ', / from eq. ( J3c| ) with t<j = 0. Points: Torque of logarithmic spirals ( Section 3.3 and Appendix C). 



sign(fc)W(q ) 



ka 



+ w(w - 20)- 



(30) 



In the second line, we have used the dispersion relation (p28|) to eliminate the derivatives of IA and W. Strictly, the term 
g/ka is negligible in the tightly- wound limit. From eqs. ( |l5| ) and (^), the fractional surface density fluctuation becomes 
<5E/E — qo — kae for go "C 1. Therefore in the linear (and stationary and tightly-wound) limit, 2iraf reduces to T WKB 
(eq. §)■ 

For stationary waves, V = 2naf is independent of radius. This quantity contains the grouping of disk parameters E 2 a 3 , 
which scales with radius as a 1-2 * 3 . Therefore if (3 = 1/2, a stationary wave occupies a single point (qo,ka) on a dispersion 
curve such as those in Fig. 1 at all radii. If j3 7^ 1/2, a radially-propagating wave train moves along the dispersion curve to 
compensate for the changes in the disk. This defines go and e as functions of a, with the angular momentum flux F as a 
parameter. In particular, if (3 < 1/2, tightly- wound trailing waves tend to become more nonlinear as they propagate inward, 
because T/(E 2 a 3 ) tends to increase with go, as is shown for a 2 = 0.2 in Figure^. Unfortunately, as Figure demonstrates, the 
pitch angle cot~ 1 (fca) also increases with go until the tight-winding approximation can no longer be trusted. This generally 
happens well before the point of streamline crossing is reached (go = 1). 

The stationary dispersion relation (^ ) can be regarded as a quadratic equation in e with coefficients in go. The two 
solutions are analogous to the short and long branches of linear theory. Since k = qo/ae, the smaller solution for e is the short 
branch. The discriminant is 



D(q ) = 



ago 



-Aga 2 q -^-{qlU) 
ago 



(31) 



The eccentricity e is real only if D(qo) > 0. Since the derivative of U diverges faster than the derivative of W as go — > 1, there 
is some go < 1 beyond which no real solutions for the eccentricity exist. At this critical go, the short and long branches join. 
These conclusions are academic, however, since the join occurs outside the tightly-wound regime. In fact, the long branch 
does not actually exist for stationary waves ( Section 3). 

We have tacitly assumed that short-branch stationary trailing waves always propagate inwards. To verify this assumption 
for our nonlinear wave trains, we must investigate the nonlinear generalization of the group velocity. An obvious candidate is 
the ratio f/p of wave flux to wave density, which reduces to the radial group velocity in the linear limit. Nonlinearly, however, 
there are actually two distinct characteristic velocities at which information about the wave propagates; neither coincides 
with f/p or Wgroup except as go — > 0. As shown by Whitham (1974), this is a generic property of nonlinear wave trains. 
The characteristic velocities in the present case are derived in Appendix B and exemplified in Figure ^. Both characteristic 
velocities are negative (inward) for trailing waves, at least as long as the wave remains tightly wound. 
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Figure 3. The two characteristic velocities (radial propagation velocities) of stationary, trailing, single-armed spirals in the tight- winding 
approximation, for f3 = 1/2, 7 = 4/3, cr 2 = 0.2 {solid curves) and <r 2 = 0.1 (dashed). The two velocities converge to the linear group 
velocity as the streamline-crossing parameter qo — > 0. Labeled points mark radial wavenumber ka along each sequence. Smaller a 1 makes 
for tighter winding (larger ka) at a given qo- 

3 EXACT SELF-SIMILAR SPIRALS 

The nonlinear wave-action formalism of Section 2 is flexible but assumes tightly- wrapped waves, \kr\ 2> 1. In this section, we 
sacrifice generality for accuracy and find solutions in the form of exact logarithmic spirals of arbitrary pitch angle. 
The surface density is 

E(r,6») =r- s/2 S{6 + alnr), (32) 

where 5* is periodic [S(ip + 2ir) = S(ip)] but may be nonsinusoidal. The cotangent of the pitch angle of the spiral, a, corresponds 
to the previous ka. As suggested by the discussion of flux conservation in section Section 2, the surface density must vary as 
r -3 / 2 on average if the spiral is to have the same strength <5E/E at all radii. Self-similarity gives the same requirements as in 
Section 2; that the sound speed vary as r _1//2 , and the adiabatic index must be 7 = 4/3. 

Syer & Tremaine (1996) have studied self-similar nonaxisymmetric power-law disks. Their solutions omit the central point 
mass and assume a constant ratio of sound speed to orbital speed (at a given phase in the spiral), as is required by strict 
self-similarity. In our case, this ratio is not constant: it varies as r-t 1- ' 3 ''' 2 in general, and as r 1//4 in this section. Nevertheless 
our models are self-similar when viewed in the limit c/v OI b — > 0, or equivalently r 2 E(r)/M — > 0, and it is only in this limit 
that they are exact. We imagine that we are studying waves very close to the central mass. 

3.1 Linear theory 

Much insight into the nonlinear dispersion relation can be obtained from the linear limit, where the surface density function 
( [32| ) reduces to 

E(r,6>) =S o r- 3/2 [l + ecos(0 + aTnr)], E (r) + eZi(r, 0), (33) 



with So a constant and e«1. The potential corresponding to (|33J) is [Kalnajs 1971 1 

$(r, 6) = -2TvGS r- 1/2 [K (0, 0) + eK(a, 1) cos((9 + a lnr)] , = $0 (r) + e$i (r, 6) , (34) 
where 



„, N 1 „ /2m + 1 + 2iq\ /2m + 1 - 2ia 
K{a,m)=-T\ K j T ( 



4 



/ 2m + 3 + 2ia\ f 2m + 3 - 2ia 



\ 4 ) \ 4 



(35) 



Linearization of the standard inviscid equations of motion yields the following dispersion relation: 
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(mil 



27rGE (r 



(r) 



«* + |Uo. 



(36) 



In the tight-winding limit a>m> 1, K(a,m) a 
pressure on the equilibrium rotation curve, 

7rG2jo(r) 2 
(cf. eq. ^|), so that the dispersion relation for m 
K(Q, 0) - 3cr 2 + (4a 2 + 9) [a 2 - K(a, 1)] = 0. 



\a\ 1 , so that eq. (|jj) is recovered. Allowing for the effects of the effects of 



1 and u) = can be put in dimensionless form: 



(37) 



(38) 



For any choice of a , one might expect two positive roots for a, corresponding to the intersection of the long and short 
branches with u = O.^j However, equation (^) has at most one root. Whereas the tight-winding equivalent eq. (^) is positive 
as kr — » 0, the left hand side of eq. ( |38[ ) remains negative as a — > because of the term involving K(a, 1). Self-gravity is 
apparently too strong to allow the long- branch wave to be stationary (see below). 

The dispersion relation has no real roots at all if a 2 J> 0.64136. At least for (3 = 1/2 and 7 = 4/3, the maximum Jean's 
length that permits these single-armed stationary waves to exist is 



(¥). 



2tt(o 



= 4.030. 



(39) 



It is clear that the real root for a, where it exists, belongs to the short branch. Since the long-branch wave depends 



primarily on self-gravity, its wavenumber should be independent of a 



0, yet eq. 



dictates that a 



a 2 — > 0. [K(a,m) f» jaj~ 1//2 for a ^> 1.] However, the wave is "short" only in a relative sense. For m 7^ 1 and except near 
resonances, the short wave normally has a wavelength comparable to the disk thickness or smaller. In this case, the typical 
wavelength is ~ A./, which is large compared to the disk thickness in the limit that Mdisk -C Af sta r (or equivalently Q — * 00) 
at fixed a 2 . For our purposes, the thin-disk approximation is almost always appropriate. 

To recap, there is no stationary long wave. However, there exist long waves with prograde pattern speeds. Equation ( |3t| ) 
shows that k 2 — D. 2 oc r -5 ^ 2 , whereas D. oc r~ 3 ^ 2 . To allow a long wave at small a, the dispersion relation (^) requires lj > 
and u) oc r" 1 . As r — *• 0, u> is asymptotically negligible compared to f2(r) but asymptotically infinite compared to the orbital 
frequency of any perturbing mass that might launch waves inward. 



3.2 Nonlinear dispersion relation 

In the limit that the self-gravity and pressure of the disk are very weak, every fluid element must follow an orbit compatible 
with the dominant potential of the central mass. Hence the streamlines must be keplerian ellipses, although their eccentricity 
need not be small. The pressure and self-gravity of the disk control the choice of the eccentricity and spirality (e, a) in a manner 
described by the dispersion relation. Since our problem is no longer local in radius, we cannot use the elegant wave-action 
methods exploited in Section 2 and instead must look elsewhere for a physical condition to provide the dispersion relation. 

Note first of all that the net torque per unit mass on each streamline must vanish, since the orbits are closed [|] But 
this condition does not provide the dispersion relation, because it can be shown that the net torques on a streamline due to 
pressure and gravity vanish separately for all choices of e and a. This is a property of logarithmic spirals for which the surface 
density varies as r~ 3 ^ 2 and the sound speed as r~ x ^ 2 along a wave crest. It is closely related to the fact that the angular 



momentum flux carried by such spirals is independent of radius [cf. syer & Tremaine (1996) and our Appendix B]. 

The actual condition from which we derive the dispersion relation is that the streamlines should not precess. Streamlines 
become free-particle trajectories if one incorporates the pressure and self-gravity of the disk into an effective potential, 

V(r,0) = -— + $(r,0) + H(r,O), (40) 
r 

where the first term on the right is due to the central mass, the second term is the gravitational potential of the disk, and the 
third term is the enthalpy, 

H =7^T> (41) 



Lyndcn-Bell & Ostriker (1967) s Antispiral Theorem requires that all roots occur in pairs ±ct: every trailing spiral (a > 0) has a 
leading counterpart (a < 0). We count only positive roots. 

t Shocks and other forms of dissipation could allow the gas to absorb negative angular momentum from the waves, and then there would 
be a nonzero torque balanced by accretion. 
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which varies with the surface density as E 7_1 . (The adiabatic exponent appears in symbolic form for the sake of clarity, but 
7 — ► 4/3.) The accelerations due to pressure are — V-ff. 

Let the time average of V along a keplerian ellipse be V(a, e, <j>), where <p is the periapse angle. Then by standard methods 
of secular perturbation theory, the condition that the orbit not precess is (cf. Brouwer & Clemence 1961) 

dV 

— (a,e,4>)=0, (42) 

to which the keplerian part does not contribute. 

It is convenient to introduce the eccentric anomaly E and the true anomaly ip. The relation between the two anomalies 
and the time in a keplerian orbit is (cf. Brouwer & Clemence 1961) 

E - esinS = ip = n(t - t pcri ). (43) 

Here n = 2n(a 3 /GM) 1 ^ 2 is the mean motion, and t per i is the time of pericenter. If <p is the azimuth of pericenter, the position 
in polar coordinates (r, 9) is 

r — a(l — e cos E), 
r cos[# — 4>{a)] = a(cosE — e), 

r sin[6> — <p{a)] = ay/l — e 2 sin£\ (44) 
Because of self-similarity, the disk potential and enthalpy take the forms 

$(a, E) = a~ 1/2 f{E), H(a, E) = a' 1/2 h{E), (45) 

where f(E) and h(E) are periodic functions. It is sufficient to calculate / and h along a single streamline. 
One can calculate f(E) as a double integral over all other streamlines: 

OO 27T 

f(E) = -a 1/2 G f da' a'T, {a') [ dE' 1 ^ eC ° sE ' 



R 



a'=0 



If 



a 2 (l — e cos E) 2 + a' 2 (l — ecos E 1 ) 2 — 2aa | (cos E' — e)(cos E — e) + (1 — e 2 ) sin E sin E' cos(a In a) 

- [{cosE' - e')sin£- (cos£- e)sinE'] yjl - e 2 sin(alna')| . (46) 

Here R is the distance between the points corresponding to E and E' on the streamlines a and a' , respectively. The factor 
(1 — ecosE') is proportional to the time spent in the interval dE' , and So(a) is a lagrangian surface density defined so that 
the mass between streamlines a and a + da as 2nT,o(a)ada — 2-7rEo(l)a -1/ ' 2 da. 

A more efficient way to calculate the gravitational potential is to decompose the surface density E(r, 6) into a sum of 
azimuthal Fourier harmonics, each itself a logarithmic spiral, and then to use Kalnajs' formula to obtain the corresponding 
harmonics of $(r, 8): 

OO OO 

E(r,0) = r~ 3/2 ^ E(m)cosm(6» + Qlnr), $(r,6>) = -27rGr~ 1/2 ^ K{ma, m)E(m) cos m(8 + a lnr). (47) 

m— m— 

Since E(r, 9) has the form (^), the angular Fourier transform needed to find the coefficients E(m) need be performed only 
at one radius. But we did not use this method except as a check. 

The enthalpy and surface density are easily calculated in streamline coordinates. The transformation (^) from (a, ip) — > 
(r, 9), when e is fixed and 4>{a) = —a In a, has the Jacobian 

_ rdrd8 r . _, ,, . 

J = , , , = VI - e 2 + aesmE. (48) 
adadip 

Since 

dip r 
— = 1 - ecosE = -, 
dE a 

we also have drd8 = JdadE. If the wave is an adiabatic deformation of a circular disk, the gas is uniformly distributed with 
respect to orbital phas ip, so that Y,rdrd8 = Y>oadadip, whence 

E(a,E) = T, (a)J-\ H(a,E) = a 2 ■ 2-KGT, {a)aJ 1 - 1 . (49) 
Now J has zeros and E has poles unless the streamline-crossing parameter 
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q = ae/\/l-e 2 (50) 

is less than unity. So q replaces the quantity qo = kae of the tightly-wrapped theory. 

The machinery above allows one to calculate the averages of <£, H, and hence V around the actual streamlines. For 
example, 

271 2n 

5>f^ = J ^(1 - e cos E)<f(a,E). (51) 
o 

But the streamline on which the average above is taken has the same shape as all of the others, whereas to compute the 
variation ( p^ ) , one needs to vary the eccentricity of the "test" streamline while holding all the rest — the "field" streamlines — 
fixed. The potentials $, H, and V are regarded as fixed functions of the Eulerian coordinates (at,Et) based on the field 
streamlines, all of which have eccentricity a and spirality a = —d<f>i/d\xiaf. the subscripts "f" and "t" indicate field and test 
quantities, respectively. A calculation yields 



(dci{\ a(e + cosE) fdEf\ _ sin E — ay/1 — e 2 cos E 

in which (at,f , -Et,f , (j>t,t ) — » (a., E, 0(a)) after the differentiation. Hence 



9$ -1/2 / dE 

-n~ — a 

oe 



J S [l (e + cos E) f{E) + ( sin E ~ a ^ 1 - e2 cos E ) je 



1 -*™ E (53) 
Jy/T^ 



The precession rate due to pressure is given by an expression identical to the above except that H & h replace <3> &/ . 

The results for a 2 = 0.2 are shown by the points in Fig. 1. We have checked the results for e C 1 against the linear theory 
of Section 3.1: a w 3.794 (linear); a w 3.792 at e = 0.01 (nonlinear). The dispersion relation predicted by the tight-winding 
theory of Section 2 is also included in Figure 1 using qo = kae rather than q as the abscissa. Both the present self-similar 
solutions and the tight-winding theory of Section 2 predict that the spirals unwind with increasing nonlinearity (da/dq < 0). 
But the long branch predicted by the tight-winding approximation does not exist. Notice that like the tight-winding dispersion 
relation, the self-similar sequence reaches a maximum q. The linear analysis of Section 3.1 shows, however, that the self-similar 
sequence cannot return to q — = e at nonzero wavenumber a. Therefore it must terminate in a point at a = 0, q > 0, and 
e — 1. The last (lowest) self-similar model shown has e = 0.998, q = 0.9016, and a — 0.057, which corresponds to a pitch 
angle of 86.7°. Our numerical methods are unable to follow the sequence beyond this point. 



3.3 Angular momentum flux 

There is more than one way to define the angular momentum flux of a spiral wave. Lynden-Bell & Kalnajs (1972, henceforth 
LBK) imagine dividing the disk in two by a cylinder coaxial with the rotation axis and consider the total torque exerted by 
the interior of the cylinder on its exterior. In collisionless disks, LBK's torque consists of two parts: a gravitational stress 
that is quadratic in gradients of the potential, and a "lorry transport" term that involves correlations between the radial and 
azimuthal velocities of stars crossing the cylinder. LBK's angular momentum flux is Eulerian — it describes the transfer of 
angular momentum from one spatial region to another. 

Lagrangian approaches define the flux by the rate of transfer of angular momentum from one mass element to another. 
Lagrangian fluxes arise naturally via Noether's theorem from an action principle in which the dynamical variables follow the 
mass, but can also be defined directly. 

For our stationary spiral waves, we define the angular momentum flux to be the total torque exerted by the disk interior 
to a given streamline on the exterior. No "lorry transport" term occurs because the gas does not cross streamlines, but the 
torque consists of two parts, a gravitational term and a pressure term. Appendix C shows that these contributions to the 
flux can be obtained conveniently by differentiating the internal and gravitational energies of the disk with respect to radial 
wavenumber aj^j Thus the prescription for the flux is the same as in the tight-winding theory; however, the gravitational flux 
so obtained is strictly accurate only for exact logarithmic spirals in disks where Eo(a) oc a _3//2 . 

For as much of the stationary sequence as we have computed, the total angular momentum flux of the trailing self-similar 
spirals is always positive and increases monotonically with streamline eccentricity (Fig. 2). 

t The total internal and gravitational energies of the disk are infinite, but we use the energies per unit logarithmic interval in semimajor 
axis, which are finite. 
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4 SUMMARY AND DISCUSSION 

We have analysed stationary single-armed spiral density waves in disks whose rotation curve is strongly dominated by the 
potential of a pointlike central mass. Both linear and nonlinear waves have been considered. For tightly-wound waves, we have 
applied a phase-averaged variational formalism to derive the nonlinear dispersion relation, angular momentum density and 
flux, and characteristic velocities of these waves ( Section 2 and Appendices A & B). The variational approach may be useful 



for other nonlinear dispersive waves of astrophysical interest (Whitham 1974). Unfortunately, when the waves are strongly 



nonlinear or when the disk is barely cool enough to permit their existence [eqs. (|2|) & (|39|)], the pitch angle is too large 
for the tightwinding approximation. Therefore we have constructed radially self-similar logarithmic spiral solutions without 
restriction on the pitch angle, but only for disks with azimuthally-averaged surface-density profiles scaling as r~ 3 ^ 2 ( Section 
3). Where both are valid, the self-similar and tight- winding theories agree well. The most important lessons learned are as 
follows (unless otherwise noted, these statements apply to stationary, single-armed waves only): 

(i) Even though self-gravity is essential, the Toomre Q parameter is irrelevant to these waves. This is because Q depends 
upon the central mass (via the orbital or epicyclic frequency), whereas the existence — and even the wavelength or pitch 
angle — of these waves is asymptotically independent of the mass ratio M centra i /A/disk, since the waves consist of streamlines 
approximating keplerian ellipses. The appropriate dimensionless ratio of temperature to self-gravity is not Q 2 but rather a 2 
[eq. (0)], which is equivalent to the ratio of Jeans length to radius. A corollary is that the waves can propagate arbitrarily 
close to the central mass provided a 2 remains below its critical value (p9[). 

(ii) These waves are nonlinear extensions of the short branch of the WKB dispersion relation (at zero pattern speed and 
azimuthal wavenumber m = 1); the long branch does not exist, at least not in isentropic disks with surface densities E oc r~ 3//2 
( Section 3.1). Consequently, trailing waves propagate inward, even nonlinearly (Appendix B). Of course the long wave does 
exist at other arm multiplicities or retrograde pattern speeds. 

(iii) If the azimuthally-averaged surface density profile is shallower (steeper) than r~ 3 ^ 2 , then ingoing (outgoing) waves 
become increasingly nonlinear as they propagate. This is a rather general result; it follows from the fact that the gravitational 
angular momentum flux scales as GS 2 r 3 at fixed pitch angle and eccentricity. 

Regarding the first point above, it is always easier to have a 2 < 0.64 in a keplerian disk than to have Q ^ 1 since 

o 2 = |-Q, (54) 

where H = v B /Q is the disk thickness. Of course the single-armed waves we have studied are not locally unstable modes. 
On the other hand, it is interesting that the Jeans length Aj = 2nra 2 is approximately the wavelength of fastest growth for 



axisymmetric secular instabilities (Lynden-Bell & Pringle 1974, Fridman & Polyachenko 1984). Viscous effects violate the 



conservation of vorticity, which tends to stabilize the disk at long wavelengths. Secular instability does not actually occur at 
large Q, however, unless the viscous torque is a constant or decreasing function of surface density ( [Schmit fc Tscharnuter 



1995) 



With regard to the second point, it should be noted that although the waves we study belong formally to the short branch, 
their wavelengths are nevertheless long compared to the disk thickness when Q S> 1. In fact their characteristic wavelength is 
Aj = ttQH. 

It is interesting that the minimum mass solar nebula, defined by augmenting the pres ent planets with enough hydrogen 



and helium to reach solar abundance, has an r 3//2 surface density profile. For example, Hayashi, Nakazawa, fc Nakagawa 



1985 quote 



E W = 1.7 X ID 3 (JL-) ~ 3/2 g cm" 3 , T(r) = 280 (^) ^ (^) ^ . (55) 

Very possibly this is only a coincidence. On the other hand, if density waves of the sort we have analysed were to dominate 
angular momentum transport, they might drive the surface density towards this power law: for a shallower profile, ingoing 
trailing waves would tend to steepen and dissipate at small radii, depositing negative angular momentum and driving accretion. 
The temperature profile above is shallower than expected for a thin optically thick disk warmed either by accretion or by 
reprocessing of solar radiation (T oc r~ 3 ^ 4 ); however, the r -1 / 2 scaling seems required to match integrated infrared spectra 



(Adams, Lada, & Shu 1987) and may result from flaring of the disk (Kenyon & Hartmann 1987; Chiang & Goldreich 1997) 



For the profiles (|55j), the parameter a is independent of radius: a w 1.5(L*/L©) ' . This is comparable to the critical value 
0.64. Under the same conditions, Q « 70(r/AU)~ 1/4 . 



VLBI observations of ex tragalactic H2O wate r masers have revealed molecular disks in the nuclei of NGC 4258 (Miyoshi 
et al. 1995| ) and NGC 1068 flGreenhill et al. 1996| ) on subparsec scales. In the former case, the disk appears to be very thin, 
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though warped, and the rotation accurately keplerian, being dominated by a 3.6 x 10 7 Mq black hole. If one assumes that 
the observed Xray luminosity of the nucleus, Lx — 4 x 10 40 erg s" 1 (Makishima et al. 1994), is powered by steady accretion 



through this disk, then (cf. Neufeld & Maloney 199E 



7 x lO-VoiaTaVai 2 . 



Q 



14.eo.oi ctT^g r o . i , 



where lOOeo.oi is the efficiency with which mass is converted to Xrays, a ^ 1 is the usual viscosity parameter, T300 is the 
midplane temperature in units of 300 K (the minimum for masing), and lOro.i pc is the distance from the black hole. Given 
the uncertainty in the first three parameters, it could easily be that Q ^ 1, but in any case the single-armed waves are surely 
permitted; the question is whether there is any way to excite them. 
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APPENDIX A: AVERAGES OVER ORBITAL PHASE FOR TIGHTLY- WOUND WAVES 

After substitution from eq. (En)), the phase average of the internal energy ( |l6|) becomes 
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2nada 



Xc 2 



7(7 - 1) J 2tt 



jT^[(l-*dn*)-^-l] =/ 



2ixada 



Ec 2 



goW(go) 



(Al) 



The phase-averaged expression is independent of 0, integration over which has yielded the factor of 2n. The definition (gj) of 
U(qo) is closely related to an integral representation of the hypergeometric function (Abramowitz & Stegun 1972), whence 



U{ql) 



J7 1 7 , 2 

aF U~2~'2' l59b 



7(7-1)^ X ' ' '"' ' ' (A2) 

In the gravitational integral (|l7|), we replace the limits on Aa = a\ — a-2 with ±00 (a good approximation for tightly- wound 
waves) and eliminate 8 in favor of *]/ using (fLSh : 



2G / aT, 2 (a)da / dAa 



d*ln 







1 + 



aecos(* + kAa/2) - aecos(* - fcAa/2) 



Aa 



The integration over ^ can be cast in the form 



27 



In II — J? sin* I 



In 



(A3) 



(A4) 



where 7? = 2ae sin(kAa) / Aa, or equivalently R = qo sin(qoy)/y if j/ = Aa/2ae. The integration over Aa can then be expressed 
in terms of the function W(go) defined in eq. (Ea). The gravitational energy is therefore 



2irada 



ae\q \W(qo) 



(A5) 



We have not been able to relate W(qo) to familiar special functions, so we evaluate it by expanding the the integrand of 
eq. (123) in powers of Oq an d integrating term by term. 



APPENDIX B: NONLINEAR CHARACTERISTIC VELOCITIES 



The characteristic velocities of a nonlinear wave train represent the speeds of propagation for modulations of the amplitude, 
wavelength, and frequency. Their derivation is simplified in our case by the special structure of the action (|2^)-(p7]) in the 
nearly stationary limit: 



I = 

Co = 



[Co{A,k,a) - LuA]J(a)dadt +0(uj 2 ), 



7rGEoa 
2aQ, 



-ge - |qo|W(go)e + a 2 qll4{qo) 



(Bl) 
(B2) 



where A = e 2 /2 and J(a) = 2-7ra 3 Eo(a)f2(a). The 0(uj 2 ) terms do not not contribute to the characteristic velocities as to — > 
and will be neglected. Variation of I with respect to A yields the dispersion relation 

dCp 
OA ' 

and variation with respect to the phase function * yields the conservation of wave action (cf. eq. 

a , T „, , d ( T ac s 



oj(A, k, a) 



(B3) 



dt {JA) + da- ( J ^k) - °- 
On the other hand, because partial derivatives of * commute, it follows from eqs. ([]) that 

dk dcu 
at aa 



(B4) 



(B5) 



The derivatives with respect to a in equations (B4) and (B5) include the indirect dependence of Co and to on a via their 
dependence on k and A. Substituting from eq. (B3) and expanding derivatives, we have 



~dt 



+ 



d 2 C \ dk 



OA 
~dt 



+ 



dAdk J da 
d 2 C \ dk 



dk 2 I da 



+ 



+ 





OA 


dCp 




da 


da 


d 2 c \ 


dA 


d 2 C 


dkdA I 


da 


dadk 



da da 



(B6) 



where now dCo/da involves only the explicit dependence of Co(A, k,a) on its third argument. 
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Equations (B6) form a system of first-order quasilinear partial differential equations for (k, A) in the (a,t) plane. Such 

V(k, A, a, t). The characteristic velocities V are 



systems can be solved by integrating along characteristic curves da/dt 
eigenvalues of the matrix of coefficients of (§f , §^-) : 



Vi,V 2 = 



d 2 c 

dkdA 



d 2 C 2 C 
OA 2 dk 2 



c)^ 
Ok. 



duj d 2 Co 
dA dk 2 



(B7) 



The reduced lagrangian (B2) has a Taylor series in A beginning at O(A), and uj has a Taylor series beginning at 0(1). 
Therefore as A — > 0, the characteristic velocities merge and can be identified with the linear group velocity, but at finite 
amplitude they are distinguished by corrections of order A 1 ^ 2 = O(e). Should the argument of the surd in (B7) be negative, 
the characteristic velocities are complex; this may signal a tendency of the nonlinear wave train to break up into a series of 
solitary waves (Whitham 1974). The required derivatives of Co are 



d 2 C 
dkdA 
d 2 C 
OA 2 
d 2 £ 



~dk 

du> 
dA 



ttGE 
tvGT, k 2 a 

In ^~ 

7rGE , 



[-a 2 ka {qlu" + 5q U' + 4W) + (q'oW + 4qoW' + 2W) sign(fe)] 



[-a 2 {qlu" + 3q W) + e (q W" + 3W') sign(fe)] 



dk 2 



29 



+ 4q U' + 2U) + e (q W" + 2W') sign(fe)] 



(B8) 



The characteristic velocities are shown for representative values of a 2 in Figure |S| We have plotted results for the short 
branch only, since the the tight-winding approximation is not applicable to the long branch. In all cases that we have examined, 
both chararacteristic velocities are real and negative for tightly wrapped trailing waves (ka S> 1, k > 0.) 



APPENDIX C: ANGULAR MOMENTUM FLUX OF SELF-SIMILAR SPIRALS 

We relate the angular momentum flux of nonlinear logarithmic spirals to derivatives of the gravitational and internal energies 
with respect to pitch angle. The flux proves to be constant in isentropic disks with Eo(a) oc a~ 3//2 and adiabatic index 7 = 4/3, 
so that no net torque is exerted on streamlines. 

Consider first the gravitational interaction energy between two keplerian elliptical rings of the same eccentricity e, 
semimajor axes (01,02), periapse angles <j>i,(j)2, and masses (mi, 7712). After integration over the arc lengths, the energy 
can be put into the form [cf. eqs (^) & (p|)] 

Wi2 = _(Gm^\ K ^ <h _ <hl (C1) 



V V" 1 " 2 / V Ol 

Because we have factored out (ai<22) -1 , the dimensionless kernel K(£,(f>) is independent of the absolute dimensions of the 
rings, and it is symmetric in the sign of both arguments. The torque exerted by ring 2 on ring 1 is, with x = In a, 

Fl2 = JT7 — = / K<t,(x2 - Xi,(f>2 - 4>i), (C2) 

o4>i v aiaa 
where K^^, <j>) = dK^, (f>)/d(j). In a continuous elliptical disk, let 

H = 27ra 3/2 E (a), (C3) 

so that the mass between streamlines a and a + da is o}l 2 \judx, and fi is constant since Eo oc a~ 3 ^ 2 . The gravitational energy 
per logarithmic interval in semimajor axis is therefore 

W{xi,<h.)=-^- I K{X2-X U 4>2-<t>l)dX2, (C4) 



2 

— OO 

and the total torque per unit x is 

oc- 

— -r(xi, 0i ) = -—jr- / Kd,(x2 — xi, tj>2 — <j>i)dx2. (C5) 
091 2 / 



In a logarithmic spiral, 4>2~(f>i = —a - (X2 — xi). Then the energy (C4) is independent of xi . However, the energy does depend 
upon the pitch angle: 
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dW - , 



The total torque exerted by all rings X2 < x on all rings Xi > x does not vanish. Referring to eq. (C2), and noting the 
antisymmetry of in its second argument, this torque is 

OO X 

r gra v(^,a) = -\~GlI 2 / dx! / dx 2 K ( f ) [xi — X 2 , — Ol{X\ — X 2 )] 



= G^ 2 J dt J driK^-at) = Gp 2 J tK*(t,-at)di = (07) 

B-£/2 

This establishes the desired relationship between the gravitational contribution to the angular momentum flux and the 
derivative of the gravitational energy with respect to pitch angle. Since dF sla , v /dx = 0, there is no net torque on streamlines. 
Now consider the pressure term. For a polytropic gas, the internal energy per unit mass is 

7-1 7-l\Vo) ' 

where p is the pressure, V = 1/E is the specific volume, and the subscript "0" denotes a reference state on the gas adiabat. 
Integrating around a streamline, and noting the area element a 2 JdEdx, we have for the internal energy per logarithmic 
interval of semimajor axis, 

2tt 

U(a,a) = ^K fj^dE, (C8) 
7-1 J 
o 

where E is the eccentric anomaly and J is the Jacobian (fes). The derivative with respect to a is 



2tt 

po(a)a 2 / esmEJ-^dE. (C9) 



dU a 
da 



The angular momentum flux due to pressure is the integral of reg T n = rpee ■ n along the arc length of the streamline, 
where T is the pressure tensor, e$ a unit azimuthal vector, and n the outward unit normal to the streamline. If a = a(r, 9) 
then n = Va| _1 Va, and the area element can be written \'Va\~ 1 dsda = JadadE, ds being the element of arc length. Hence, 

r prcss = / P {a,E)aJ dE. (CIO) 



Now jp(a, E) — T,o(a)c 2 i o(a)J^ 7 , and equations ( pfi| ) yield (da/dd) r = — ^7 1 aesin£'. With these substitutions, we have 

Tpress — T: , (Cll) 

oa. 

as claimed. Since the prefactor p^{a)a 2 is independent of a in our self-similar spiral disks, r pross is also independent of a, so 
that the net pressure torque on a streamline vanishes. 

This paper has been produced using the Royal Astronomical Society/Blackwell Science F/IpjX style file. 
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